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ABSTRACT 



We report the detection in Ks-band of the secondary echpse of the hot 
Jupiter CoRoT-lb, from time series photometry with the ARC 3.5-m telescope 
at Apache Point Observatory. The echpse shows a depth of 0.336±0.042 percent 
and is centered at phase 0.50221q;qq27; consistent with a zero eccentricity orbit (e 
costu = 0.0035^q;qq42). We perform the first optical to near-infrared multi-band 
photometric analysis of an exoplanet's atmosphere and constrain the reflected and 
thermal emissions by combining our result with the re c ent 0.6, 0.71, a nd 2. 09 /im 



secondary ec li pse de tections by iSnellen et al 



Alonso et al. 



fl2009h 



Gillon et al. 



mm . and 



(j2009al ). Comparing the multi- wavelength detections to state-of- 
the-art radiative-convective chemical-equilibrium atmosphere models, we find 
the near-infrared fluxes difficult to reproduce. The closest blackbody-based and 
physical models provide the following atmosphere parameters: a temperature 
T = 2454+^^0 K, a very low Bond albedo Ab = 0.000+°;°^^, and an energy re- 
distribution parameter P„ = 0.1, indicating a small but nonzero amount of heat 
transfer from the day- to night-side. The best physical model suggests a ther- 
mal inversion layer with an extra optical absorber of opacity Kg = 0.05 cm^ g~^, 
placed near the 0.1-bar atmospheric pressure level. This inversion layer is lo- 
cated ten times deeper in the atmosphere than the absorbers used in models to 
fit mid-infrared Spitzer detections of other irradiated hot Jupiters. 



Subject headings: binaries:eclipsing - planetary systems - stars:individual (CoRoT-1) 
- techniques: photometric 
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Introduction 



Space-based detections of liot Jupiter atmospheres have flourished in recent years. The 



at wavelengths longer than 3.6 fj,m (e.g. 


Charbonneau et al. 


2005 




Deming et al. 


2005 




2007; 


Knutson et al. 


2007. 


2008; 


Harrington et al. 


2007; 


Machalek et al. 


2008. 


200 


9), while 



near-infrared (1.5-2.5 /im) observations with the Hubble Space Telescope have found 
evidence for water, carb on monoxide, and carbon dioxide in the dayside spectrum of the 



exoplanet HD 189733b (ISwain et al. 



20091 ). These detections have been made during 
secondary eclipses (when the planets pass behind their host stars). Important atmospheric 
absorption signatures have been detected by transit observa tions with these space telescopes 



as well: sodium in HD 209458b (ICharbonneau et al. 



(ISwain et al. 



20071 1. 



2OO2I) and methane in HP 189733b 



20081 ) with Hubble, and water vapor in HD 189733b with Spitzer (ITinetti et al. 



The CoRoT mission recently joined these space-borne successes by detecting phase 
brightness variations and comb ined thermal and reflected emission duri ng secondary eclipses 



of the exoplanets Co RoT-lb (ISnellen et al. 



2009 



Alonso et al. 



( lAlonso et al. 



2009aJ) and CoRoT-2b 



2009bl ) in an optical broadband window centered at about 0.6 fim (see 



Table [I]). 

At the same time that these detections have provided very valuable insights into the 
atmospheric physics of irradiated hot Jupiters, they have also revealed some perplexing 
flndings. Some of t he planets show a strong temperature cont rast between their day 



and night sides (e.g. 



Harrington et al 



2006 



Snellen et al. 



have a more efficient redistribution of incident energy (e.g. 



20091). while others a ppear to 



Knutson et al. 



20071 ). Another 



striking discovery is an apparent bifurcation of hot J upiters into two clas s es based on the 



presence or absence of a thermal inversion layer (e.g. 



Hubeny et al. 



2003 



Burrows et al 
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2008bl : iFortnev et al 



[see e.g. 



Barge et al 



20081) . Along with a wider-than-expected range of exoplanet radii 



20081 ) ■ these are currently the key unsolved questions in the study of 



irradiated hot Jupiter atmospheres. 

Resolving these questions will require observations from the optical to the infrared 
wavelength regimes. Ground-based observations have just recently started to reach 
the sensitivity necessary to directly detect hot Jupiters. The first two detections were 
announc ed in early 2009 at optical a nd near-infrared wavelen gths: OGLE-TR-56b in 



z'-band (iSing fc Lopez-Morales 



20091), and TrES-3 in K-band (jde Mooij fc SnellenI 120091 ). 



he third ground- based result was the detection of an eclipse of CoRoT-lb at 2.09 /im 



( iGillon et al. 



20091 ). This last detection makes CoRoT-lb the first exoplanet to have 



thermal emission measured at both optical and near-infrared wavelengths. 

Here we report the fourth ground-based detection of thermal emission from an 



exoplanet, and the fourth detection of CoRoT-lb, this time in Ks-band (2. 
combi n e our results with t he oth er three recen t detec tions of CoRoT-lb by 



(120091) 



Gillon et al 



3) 



5 /xm). We 



Snellen et al. 



((20091), and 



Alonso et al 



(l2009al ) to produce the first multi-color 



analysis of the atmospheric spectrum of an exoplanet between 0.5 and 2.2 /xm. 

Section [2] describes the observations. In section [3] we detail our reduction and analysis 
steps to detect the signal from the planet. Section H] compares the detected signals to the 
predictions made by state-of-the-art planetary atmosphere models. The results are then 
discussed and summarized in sections [5] and El 



2. Observations 



We observed two secondary eclipse events of CoRoT-lb on the nights of 2009 January 9 
and 15, UT, using the Near-Infrared Camera & Fabry-Perot Spectrometer (NICFPS) on the 
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ARC 3.5-m telescope at Apache Point Observatory in New Mexico. Windy conditions and 
thin and constantly varying cirrus layers hampered our attempts to reach high-precision 
photometry in the January 9 observations. The conditions on January 15 produced good 
quahty data, which we will discuss here. 

NICFPS is equipped with a Rockwell Hawaii 1-RG Ikxlk HgCdTe detector with a 
4.58x4.58 arc minute field of view and a pixel scale of 0.273 arcsec/pix. We chose the 
reddest available broadband filter (Ks), in order to maximize the eclipse depth and collect 
a large number of photons in as short an exposure time as possible. The instrument has a 
high read noise (95 e~/pixel), which can be reduced through the use of up to 15 Fowler 
samples. In this mode, the chip is read out non-destructively a number of times during the 
exposure, and the difference of each pair of readouts is taken as one sample. 

We used the instrument in the standard 8 Fowler sample mode, in which the chip is 
read out 16 times consecutively while exposing. The first readout is subtracted from the 
ninth, the second from the tenth, and so forth, producing eight 5.44-second samples in 10.88 
seconds. We then use the average of the eight resulting differences, which reduces the read 
noise of the detector to 35 e^ /pixel (a factor of y/S). To achieve a good sky subtraction, we 
used a simple two-point dither pattern, taking a set of two exposures and then offsetting 
the objects by 58 pixels (15.8 arcsec). This resulted in two exposures every ~50 seconds, 
when combining the exposure, readout, and offset times. With a measured gain of 4.77 
e~/ADU, we collected ~782,000 photons per datapoint from the target, and 746,000 to 
1,417,000 photons pe r datapoint fro m the nearby comparison stars. Based on the noise 



expression derived by 



Howelll (119891 ) to combine the Poisson noise with the sky level, dark 



current, and read noise, we expected the noise per datapoint to be 0.259%. 

In order to capture the full secondary eclipse and a good baseline for the light curve, 
we monitored the target from 04:20 to 09:40 UT (twice the eclipse's duration), obtaining 



- 6 - 



758 frames over the course of 5.3 hours. While the transparency was excellent, the seeing 
varied between 0.97 to 1.91 arcsec. Throughout the observations, slow drifts in telecope 
pointing led to a varying target location within a ~4-pixel circle. 

In addition to the target observations, we took 120 dome fiats with dim quartz lamps 
in the Ks filter, and 400 darlc frames, each for a range of exposure times from 1 to 8 seconds 



3. Reduction and Analysis 

We began the data reduction by subtracting from each dome fiat the average of the 
darli frames for the fiat's exposure time, and then combining all the fiats to create one 
master normalized fiat. We then applied the fiat-field correction to the target images, 
and created reduced, sk;y-subtracted images by subtracting from each image the nearest 
neighbor frame talcen at the opposite dither position. 

Two team members then performed independent analyses of the data, as discussed in 
sections 13.11 and 13. 2[ in order to minimize systematics and achieve the optimal photometry. 
In both analyses we converted the JD in the headers of the images to HJD a nd then to a n 



orbital phase using the most recent ephemerides for transits of CoRoT-lb by 



Bean 



(120091 ). 



3.1. Analysis A 

The first approach analyzed the combined dither positions as a single dataset. Several 
stars, including the target, were isolated enough in the frames to be analyzed by standard 
aperture photometry. 

Using the IDL adaptation of DAOPHOT, we recorded the fiux from CoRoT-1, plus 
another 20 bright, isolated stars in the field, for aperture sizes from 1 to 20 pixels, in 
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increments of 0.05 pixels. Of the 20 field stars, we selected as comparisons the four that 
produced the most stable flux ratios with respect to the target, i.e. the differential light 
curves with the least photometric dispersion during the predicted out-of-eclipse phases. 
The locations of these four best comparisons are shown along with the target in Figure [H 
A fixed aperture size of 6.1 pixels (1.67 arcsec) and a 14-24 pixel (3.82-6.55 arcsec) sky 
annulus produced the most stable differential photometry when combining the curves from 
all four comparisons. The average photometric dispersion in the out-of-eclipse portion of 
the combined light curve was 0.781%, a factor of 3.02 larger than the expected noise per 
datapoint. 

However, the combined curve also showed clear systematic trends that could be 
attributed to atmospheric effects such as seeing and airmass, as well as instrumental effects 
such as changes in the (x,y) position of the stars in the images or temperature and pressure 
changes in the instrument throughout the duration of the observations. 

These trends were individually investigated by fitting a linear correlation between each 
parameter and the differential flux, for the out-of-eclipse points only. The most significant 
systematic trends correlated with variations in the seeing, and we found that nearly all 
the systematic noise in the light curves was removed after removing a trend based on that 
parameter. 

For the final photometry we de-trended each differential light curve with respect to 
seeing, and then combined the four, ending up with a dispersion of 0.547%, 2.11 times the 
expected noise limit. The final light curve from this analysis is shown in the top panel of 
Figure [2l 
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3.2. Analysis B 



The second analysis approach, analogous to that of ISing &: Lopez- Morales! ( 120091 ). 
began with the separation of the tw o dither position se ts into different light curves, and 



implemented the SysRem algorithm (ITamuz et al 



20051 ) for de-correlation. 



We used standard IDL procedures to perform aperture photometry on 19 stars, 
including the target, for aperture sizes between 1.0 and 14.9 pixels spaced in increments of 
0.1 pixel. The centers of each star were determined by fitting a 2D moffat function to the 
point spread function (PSF) of each star. The residual sky background was determined by 
finding the sky annulus which resulted in the highest target signal-to-noise (photometry 
error with photon, sky, and readnoise). On the basis of this evaluation, the best sky annulus 
resulted to be one of inner radius 19 pixels and outer radius 20 pixels, corresponding to 120 
pixels between 5.19 and 5.46 arcsec from the star. 

The five stars selected as comparisons in this analysis are shown in Figure [H The 
comparison selection criterium was the same as in Analysis A, i.e. the field stars that 
produced the most stable photometry with respect to the target during the out-of-eclipse 
phase. Later comparison showed that the four reference stars used in Analysis A overlap 
with these five. 

After producing individual light curves for each dithe r position, we de-co rrelated the 



curves using an implementation of the SysRem algorithm (iTamuz et al 
to minimize the expression 



20051 ). which seeks 



^{rij -aajf/a^, (1) 

where Vij is the average-subtracted stellar magnitude for the ith star of the jth image, a is 
the uncertainty of epoch-dependent parameter, and aj is a stellar-dependent 
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parameter. The first pass through the SysRem algorithm produced a hght curve with a 
single clear linear trend. We found this linear trend to be efficiently removed by either 
a second pass through the SysRem algorithm or by allowing the baseline ffux to vary in 
time linearly, fit by two parameters. Removing this linear trend effectively removed any 
detectable systematic trends. We fit each dither position individually at first, but found 
that the linear slope was similar for both, and re-combined the sets. 



In order to test the effectiveness of our de-correlation procedures, we searchec 



presence of residual systematic errors correlated in time ( "red-noise" , iPont et al 



for the 



20061) 



by checking that the binned residuals followed an N~^/'^ relation, when binning in time 
by points. The presence of red-noise causes the variance to follow a cr^ = cr^/A^ + 
relation, where is the uncorrelated white noise component while cr^ characterizes the 
red-noise. We found no significant evidence for red-noise when binning on time scales up to 
42 minutes (100 points), showing that any correlated noise had been effectively removed by 
the de-trending. 

The optimal aperture size of 5.6 pixels (1.53 arcsec) was determined to be that which 
minimized the standard deviation of the light curve when using the target and reference 
stars in these de-trending procedures. The final light curve from this analysis had a 
dispersion of 0.661%, 2.55 times the expected noise limit. The final light curve from this 
analysis is shown in the second panel of Figure [2j 



3.3. Secondary Eclipse Fit 



The next step in our analysis was to search for the eclipse signal from the planet in the 
observed light curves. This was done separately for each light curve resulting from analyses 
A and B, and the solutions were then combined to produce a final result. 
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Gillon et al 



We used the o rbital period from iBearu (120091 ) and the stellar and planetary radius from 



(120091 ) (see Table [3]) to generate a grid of eclipse models with no limb darkening 
and a constant out-of-eclipse baseline. The free parameters in the model grid were the 
baseline level and the depth and central phase of the eclipse, with the grid covering baseline 
levels between 99.6 and 100.4% in increments of 0.002%, eclipse depths between and 0.8% 
in 0.002% steps, and central phases between 0.46 and 0.54 in 0.0001 phase increment steps. 
Running each light curve through the model grid resulted in the following best-fit solutions: 
For the analysis A light curve, a baseline of 100.014% and an eclipse of depth 0.322%, 
centered at phase 0.5003; for the analysis B light curve, a baseline of 100.006%, an eclipse 
depth of 0.355%, and a central phase of 0.5073. We decided to use the point-by-point 
average of the light curves from analyses A and B for our final result, shown in the third 
panel of Figured Table [2] shows the first few unbinned points in the final result, and the 
full dataset is available in the online version. The best-fitting model in that case, found via 
chi-square minimization and with a reduced of 1.051, gave a baseline of lOOMOtoZl^ 
an eclipse depth of 0.336^Q'gg4 and a central phase 0.5022+f^:|][{2?- 

These values are given, 

along with the other parameters of the system, in Table El The results of each fit for the 
entire parameter space are shown in Figure [31 

To estimate the error in the eclipse depth, we put all the in-eclipse points from the 
de-trended photometry into a single bin, the entire out-of-eclipse portion into a second bin, 
and combined the binned errors. We had 295 in-eclipse points with an average individual 
dispersion of 0.531%, and 344 out-of-eclipse points with an average individual dispersion of 
0.528%, which combined to give us a la error of 0.042%. For the error in phase we have 
directly adopted the la confidence values from the contour plot in Figure [31 

As an additional test to confirm the depth of the eclipse, we generated histograms of 
the distribution of normalized flux for both the in-eclipse and out-of-eclipse portions of the 
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light curve, adopting as central eclipse phase the value ip = 0.5022 obtained above. The 
result, shown in Figure HI shows a clear 0.336% shift of the distribution of in-eclipse points 
with respect to the out-of-eclipse points, in full agreement with the result of the model grid 
fits. 



4. Atmosphere Model Fits 



While secondary e clipse detections in a single broac 



with the CoRoT data ([Snellen et al. 



2009 



Alonso et al. 



band, such as the optical detections 



2009bl ). cannot constrain a planet's 



thermal and reflected light independently of one another, we can begin to disentangle the 
reflected and thermal contributions to the total light from the planet by comparing the 
secondary eclipse depths at different wavelengths. Combining our Ks- b and eclipse detectio n 



of C oRoT-lb with the o ther three detections by iSnellen et al.l (|2009|), iGillon et al. 



and 



Alonso et al 



fl2009h . 



(I2009al ). in this section we make the first simultaneous multi-wavelength 
comparison of observations of an exoplanet atmosphere to current models at optical and 
near-infrared wavelengths. 

In the following subsections we compare the observations first to simple blackbody 
models and then to more sophisticated radiative-convective models of irradiated 
planetary atmospheres in chemical equilibrium. The advantage of the blackbody models 
is that they provide simple initial estimates ofthe global properties of the planetary 
atmospheres, such as their approximate temperatures, their reflective properties, and 
how effi ciently they redistribute energy from the irradiated to the non-irradiated sides 



[see e.g. 



Lopez-Morales fc Seager 



20071 ). Atmosphere models provide more details about 



specific properties of the atmospheres, such as finer spectral information. This allows 



the identification of specific absorbing or emitting chemicals, as well as t 
atmospheric depths, and the presence of thermal inversion layers (see e.g. 



leir specific 



Burrows et al 



2008aD. 



4.1. From Contrasts to Flux Densities 

Before beginning the model comparisons, we converted tlie measured eclipse depths 
to planetary fluxes, because secondary eclipse photometry does not provide an absolute 
measure of the planet's brightness, but a planet-to-star flux ratio integrated through the 
instrument and filter profile. 

For this conversion, we calculated the synthetic flux density of a Kurucz (1993) 
GOV-type star, equivalent to the CoRoT-lb host star, in the CoRoT White, CoRoT Red, 
NB2090, and Ks filters using each filter's response curve. The CoRoT White, NB2090, 
and Ks filters have well-defined transmission functions centered at 0.6, 2.09, and 2.15 fim, 
respectively, but the response curve for the CoRoT Red channel is not as well-known. The 
collected light is passed through a prism that divides it into red, gr e en, an d blue portions. 



Snellen et al. 



(120091 ) ■ we estimated 



but the behavior is different for each star. Following 
a wavelength cutoff of ~560 nm by taking the fraction of light in the red channel to the 
total light, giving an effective wavelength of 710 nm. To approximate the behavior of the 
red channel, we used the same response curve as for the white channel, but with zero 
transmission below 560 nm. 

For the purposes of the comparison we used a hypothetical distance of 10 pc (i.e. the 
physical flux density of the star and planet were derived for this standard scale distance, 
not their actual distance). We then multiplied each synthetic stellar flux density by the 
observed planet-to-star contrast in that band. The resulting planet flux densities are 
summarized in Table [H together with the observed contrasts and the parameters of each 
fllter. The derived planet flux densities are also shown in Figure [51 
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We note that the contrasts observed and the flux densities calculated for CoRoT-lb in 



the NB2090 and the 



^s filte rs are very similar, our measurement confirming the detection 



made by iGillon et al 



fl2009h . 



4.2. Blackbody-based Planet Models 

In the blackbody approximation, the temperature profile of an irradiated planet will be 
a smooth combination of refiected and thermally emitted light. That temperature profile is 
determined by the interplay of the stellar irradiation. Bond albedo As, and a re-radiation 
factor /, which describes how efficiently energy from incident radiation is transport ed 



Lopez-Morales fc Seager (120071), the 



3, 



around the planet before being re-emitted. Following 
value /=2/3 corresponds to no redistribution before re-radiation and /=l/4 corresponds 
to the incident energy being evenly redistributed around the planet. The combination of 
the optical and near-infrared measurements of CoRoT-lb provide powerful constraints on 
both the refiected hght (i.e. planetary albedo) and the heat redistribution between the day 
and night sides of the planet. In this section we apply models combining refiected light 
and thermal blackbody radiation to interpret the detections. The results are illustrated in 
Figures O and O 

Figure [5] includes the four measured planetary fiux densities, normalized to a distance 
of 10 pc, along with models with different temperature, albedo and re-radiation factor 
values. The effect of increasing the albedo is to increase the flux density measured in the 
optical (which corresponds to reflected light), while decreasing the infrared emission due 
to the lower absorbed energy. Increasing the day-night side heat redistribution lowers the 
entire spectral energy distribution, but does not shift the balance between the optical and 
near-infrared flux densities. The models shown in the flgure reveal that the observations 
are entirely inconsistent with efficient day-night heat redistribution and also exclude a high 
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albedo. 

Figure E] shows the result of exploring the entire - / parameter space using the 
observed flux densities and the thermal emission plus reflected light models. To quantify 
the constraints and identify the best fit, we evaluated these models along a grid spanned 
by Bond albedos Ab=0 to 0.3 and re-radiation factors f=l/4 to 2/3, calculating at each 
grid point the predicted temperature, the predicted flux densities in the four observed 
bands, and a value. The contour levels in the figure show the total four-dimensional 
distance of each model from the observations, expressed in units of uncertainties. This was 
calculated by summing the squares of the differences between the observed and predicted 
flux densities divided by the uncertainty at each wavelength, and taking the sum's square 
root, i.e. equivalent to the square root of the value. The best fit corresponds to a model 
with dayside temperature 24541;^7q K, a very low albedo {AB=OmOto fol), and inefficient 
but measurable heat redistribution (/=0.450lo Qgg). These parameters are given in the 
third section of Table [31 Thus, our simple blackbody-based modeling suggests a low-albedo 
planet with some measurable levels of energy re-distribution, but still a very prominent 
temperature difference between the day and night sides. 

4.3. Theoretical Atmosphere Models 

Due to prominent molecular absorption bands heavily-irradiated giant planet 
atmospheres are thought to display an emission spectrum very unlike a blackbody. Thus, 
while the blackbody models provide useful first estimates of the planetary properties, 
realistic atmospheric models are required to derive the actual physical properties of the 
planets. To gain more detailed information about the physical and chemical processes 
undergoing in the atmosphere of CoRoT-lb, we compared the observations to the latest, 
and still evolving, models of irradiated hot Jupiter atmospheres. 
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The model atmospheres we used are derived from self-consistent coupled radiative 



transfer and chemica l equi li brium calculat i ons, 



Sudarsky et al. 



Dased on the models described in 



(12000 



20031 ) 



Hubeny et al. 



fl2003h . and 



Burrows et al 



(12005 



2006 



2008ah . 



The most important components of the code include molecular and atomic opacities, 
and calculations to determin e the chemical abundances using thermochemical models 



[e.g. 



Sharp fc Burrows 



20071 ). The day and night sides of the planet are trea ted separately , 



with the day side receiving incident flux from the star using the appropriate 



Kuruczl fll993h 



spectral model, and the night side receiving heat from the day side via con vection. The 



convec tion is modeled with a mixing length equal to the pressure scale height (IBurrows et al. 



2008al ). The heat redistribution, desc ribed by a parameter P „,, is used to derive the planet's 



flux at the time of secondary eclipse (IBurrows et al. 



2008bl ). The P„ parameter represents 



the fraction of the incident stellar energy that is redistributed to the night side (P„=0 
corresponds to no redistribution; P„=0.5 is uniform distribution around the planet). Q 

The immediate conclusion from applying the base models described above is that 
atmospheres with no thermal inversion, unable to reproduce the near-infrared fluxes we 
observe, can be confidently excluded. Therefore we modified those models by adding an 
extra optical (0.37 to 1.0 /im) absorber with a constant opacity Kg to the abundances 
predicted by the chemical equilibrium calculations. This absorber, placed at an atmospheric 
height (i.e. pressure) of P = 0.01 bar, has the effect of creating a strong temperature 
inversion, with the extra optical absorption heating the stratosphere. Some of the resulting 
representative models are shown in Figure [71 



^ Pri=0 (no energy redistribution) corresponds to /=2/3 and P„=0.5 (maximum redis- 
tribution) corresponds to /=l/4. However, as the physical models incorporate factors that 
are accounted for differently than in the blackbody models (e.g. pressure, opacity), the Pn-f 
relation is degenerate, i.e. multiple P„ values may correspond to the same / value. 
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The upper panel in the figure shows identically irradiated planets, but with different 
redistribution parameters (P^ = 0.1, 0.3, 0.5) and extra absorber opacities [k^ = 0.0, 
0.1, 0.0 cm^ g^^)- Models with no energy redistribution (i.e. P„=0.0) were too bright in 
the optical regime to fit the two CoRoT points. Although representing diverse planetary 
atmospheres, none of these three, nor any of the other models, reproduced the observed 
very bright near-infrared fiux densities. 

In an attempt to more closely reproduce the high observed near-infrared fiuxes, we 
used an atmosphere model with P„=0.1 and Kg = 0.05 cm^ g~^. The extra absorber had to 
be placed deeper in the atmosphere (at ~0.1 bar). This is a factor of 10-100 deeper than 
what has been used to fit the Spitzer data for other hot j up iters, in whi ch case an extra 



optic al absorber at ~10 - 10 " bar reproduced the IRAC points (e.g. 



Machalek et al. 



20081 ). The temperature at this layer, also calculated by the model, is around 2200 K. 
The result of this final model attempt is illustrated in the lower panel of Figure [3 This 
model fits the NB2090 and CoRoT channel observations fairly well (within 0.9 cr), but still 
underpredicts the observed Ks-band fiux density by 2.6 cr. Still, there is some improvement 
over the upper-panel models, for which the Ks-band deviations are between 3.3 and 4.9 a. 
The parameters in this best-fit model are also listed in Table |3l 

Previous to this work, there were virtually no observational constraints on hot Jupiters 
in this wavelength range. Although a perfect model match to the near-infrared data has 
yet to be achieved, the results reveal very hot temperatures at low-to-moderate optical 
de pths, i.e. in the upper stratospheres probed by the 2-micron observations (e.g. Figure 1 



m 



Burrows et al. 



2008ah . 
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5. Discussion 



The detections of exoplanets via transits and secondary eclipses have produced many 
of the most valuable insights to their physical properties, but also yielded some puzzling 
results. Among these is the apparent division of the population of hot Jupiters into two 
distinct classes based on their atmospheric properties: one group that appears cooler, with 
water and methane absorption bands and more efficient energy redistribution, and another 
with higher levels of thermal emission, strong day-night contrasts, and lacking the expected 
absorption bands. The dominant explanation for this dichotomy is the presence in many 
planets of a thermal inversion la.yer, i.e. a hot stratosphere caused by extra absorbers of 



optical light (IHubeny et al. 



2003; 



Burrows et al. 



2008bl : 



Fortney et al. 



20081). 



The exact natu re of the absorbers is a difficult question still being investigated. 



Spiegel et al. 



20091 ) recently argued that vanadium oxide is not likely to fulfill this role, 
and that the previously favored titanium oxide would require unusually high levels of 
macroscopic mixing to remain in the upper atmosphere. S2, S3, and HS compounds, as 
absorbers of optical and ultraviolet light, hav e also recently been both considered and 

1). Co n sistent modeling 



2009a 



Burrows et al. 



fl2008ah suggest 



questioned as causes of the thermal inversion (iZahnle et a 
of the infrared secondary eclipse spectra of six planets by 
that the presence of the necessary absorber may depend not only on the incident stellar 
radiation, but also on planetary metallicity and surface gravity. 

From the combination of the four optical to near-infrared detections we have detailed 
in this paper, it appears that CoRoT-lb clearly falls into the class of hot Jupiters with a 
thermal inversion. The best-fit physical model to the combined dataset requires a small 
redistribution parameter P„ = 0.1 and an extra optical absorber with fiat opacity k = 
0.05 cm^/g. While we cannot determine the identity of this absorber, the model constrains 
its altitude, requiring absorption a factor of 10-100 times deeper in the atmosphere than 
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suggested by previous results for other exoplanets. This new constraint also highlights the 
power of combined optical and near-infrared photometry. 

CoRoT -lb reflects some other surprising trends as well. Hot Jupiters tend to have very 



low albedos ( Sudarsky et al. 



2000 



Rowe et al. 



2006 



Burrows et al. 



2008bl ) , and our findings 



place this planet's atmosphere in agreement with that trend. From the best-fit physical 



model derived in Section H73| we obtained an estimate of the planet's geometric albedo at 
the wavelengths probed by the CoRoT detections (0.4 to l.O/xm) to be Ap = 0.05±0.0. 



Assuming a wavelength-independent Lambert sphere (see 
this corresponds to a Bond albedo Ab=0.075±0.015. 



Lopez-Morales fc Seager 



2003), 



CoRoT-lb is also notable for its extremely large radius (1.45 R r,,^. iGillon et al.ll2009l ): 
evolutionary models f or a hot Jupiter of it s mass, age, and irradiation predict a radius of 



only 0.94 - 1.18 Rjup flFortnev et al. 



20071 ). This places it towards the upper end of the 



wide distribution of radii that has been seen among the population of transiting exoplanets. 

While small-radius planets can be modeled with a larger, denser core, the inflated 
sizes of planets like CoRoT-lb provide quite a challenge to planetary models. Several 
explanati ons have been propo sed, including heat retained by enhanced atmospheric 



opacities (IBurrows et al. 



( IShowman et al. 



20071 ). deposition of kinetic wind energy in the upper atmosphere 



Ibgui fc Burrows 



and significant tida 


heating caused 


orbital eccentricity or 


1 

IBodenheimer et al. 


2003 




Greenberg 


2009; 


Marlev et al. 


2009 



20091). Though 



Gillon et al 



(I2OO9I ) measured an eccentricity e=0.07llo; 



+0.042 
028 



for CoRoT-lb's orbit, the eccentricity we measure (e cosco'=0.00351q;qq42) based on its 
mid-eclipse phase is consistent within its errors with a circular orbit. Thus, there is 
no evidence that tidal heating derived from a currently eccentric orbit contributes to 
the energy bu dget of CoRoT-lb , although intense heating in the recent past cannot be 



excluded (e.g. 



Miller et al 



(I2OO9I )). In each of the diverse models, the radius of the planet 
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is determined by the pressure-temperature structure of the atmosphere and the balance 
between the energy input and output. In order to understand and improve the models, we 
must first constrain the energy balance by understanding the basic atmospheric properties, 
such as chemical abundances, pressure-temperature profiles, and energy redistribution. 

With our detection in Ks joining the 0.6, 0.71 and 2.09 /^m measurements, we are 
beginning to get a more detailed picture of the atmospheric properties of CoRoT-lb. 
However, there are still sizeable gaps in the observed spectrum, and even the most advanced 
current models have difficulty explaining the high flux levels in the 2-fim window. It is 
crucial to add further detections in both narrow and broad bands, in the optical and 
near-infrared, and a larger target sample to provide stronger constraints to the planetary 
atmosphere models. 

6. Summary 

- We directly detect thermal emission from CoRoT-lb in Ks-band, determining the 
planet-to-star flux ratio to be 0.336±0.042%. 

- Using simple blackbody-based models, we find the best fit to be a blackbody of 
2454l^yQ K, and confidently rule out both a high albedo and efficient day-night heat 
redistribution. 

- Using realistic atmosphere models, we find the need for a thermal inversion layer 
and an extra absorber near the 0.1-bar level, deeper in the atmosphere than the Spitzer 
mid-infrared data suggested for other hot jupiters. 

- Both the blackbody models and physical atmospheric models agree on a small but 
non-zero amount of heat redistribution, and a Bond albedo less than 0.09. 
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In short, the combined optical to near-infrared photometry of CoRoT-lb has allowed 
us to independently constrain the reflected light and thermal emission, and has revealed a 
very hot, low-albedo planet with a large day-night contrast and a prominent temperature 
inversion. 
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Filter 


Ae// (/iin) 


FWHM (/iin) 


Peak Trans. 


Planet: Star Flux Ratio 


PI. Flux Density (Jy 


CoRoT white 


0.60 


0.42 


72% 


0.016± 0.006%! 


0.0093±0.0035 


CoRoT red 


0.71 


0.25 


72% 


0.0126± 0.0033%2 


0.0087±0.0023 


NB 2090 


2.095 


0.020 


82% 


0.278t°;°i%3 


U.iUy^_Q 0260 


Ks 


2.147 


0.318 


97.5% 


0.336±0.042%^ 


0.1172±0.0160 



Table 1: Transmission information of each detection filter, along with the planet-to-star flux 
ratio and physical flux density of the planet in each fllter. The p hysical flux densities are 



calculated for a hy pot hetical distance o 



Snellen et al 



20091 . (3) 



Gillon et al. 



10 pc. References. (1) 



Alonso et al. 



2009al . (2) 



20091 . (4) this work. 
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Fig. 1. — Finder chart for CoRoT-1 target (enclosed in square) and the field stars selected 
as photometric comparisons. The four stars used in Analysis A are enclosed in the larger 
circles, while Analysis B used those four plus the star in the smaller circle. 
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Fig. 2. — Top Panel: Light curve from Analysis A after de-trending, with its best-fit modeL 
Second Panel: Light curve and best-fit model from Analysis B after de-trending and removal 
of the remaining linear trend. 

Third Panel: Final light curve calculated by combining Analyses A and B. A model with 
the best-fit secondary eclipse is shown as the horizontal red line, with best-fit central phase 
shift of 0.5022. The baseline is set at 100.010%, with an eclipse depth of 0.336%. 
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Center phase 
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Fig. 3. — Top Panel: Model eclipse depth versus central phase for phases between phases 
0.46 and 0.54. The best model fit to the data has a depth of 0.336% and central phase 0.5022. 
The eclipse depth falls rapidly in both directions away from that phase value. Bottom Panel: 
Confidence contours of the best fit at the 68.3%, 95.5% and 99.7% level. The best-fit value 
is indicated as a cross at phase 0.5022. The star at 0.5 indicates the phase of the expected 
center of the eclipse for a circular orbit. The result is therefore consistent with a zero 
eccentricity orbit (e cosu= 0.0035 laooD- 
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Fig. 4. — Normalized flux histograms of the in-eclipse (red dotted line) and out-of-eclipse 
(black solid line) portions of the CoRoT-lb light curve in Figure [2l The width of each bin 
is 0.336%, the same as the detected eclipse depth. 
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Fig. 5. — Comparison of the four detected CoRoT-lb planet fluxes with models that include 
reflected light and heat redistribution, for a range of Bond albedo (A^) and re-radiation 
factors (/). We show two models with maximum energy redistribution (/ = 1/4) and Bond 
albedos of (red) and 0.3 (blue), and two models with no energy redistribution (/ = 2/3) 
and Bond albedos of (green) and 0.3 (purple). The best-fit model, shown in black, is for 
a zero-albedo planet with a small, but non-zero amount of heat redistribution. 
At the bottom of the figure we show scaled transmission functions for each of the filters used: 
CoRoT the dashed line shows the blue cutoff for the red channel), NB2090 (solid), and Ks 
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Albedo and Re-radiation Model Fits 
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Fig. 6. — Contour plot showing the best-fit albedo and re-radiation factors for the planet 
fluxes detected. The contour levels show the total four-dimensional distance of each model 
from the observations, expressed in units of uncertainties. This was calculated by summing 
the squares of the differences between the observed and predicted flux densities divided 
by the uncertainty at each wavelength, and taking the sum's square root, i.e. equivalent 
to the square root of the value. The data clearly favor a very low-albedo planet with 
inefRcient but measurable energy redistribution, with the best fit (indicated by a bold 
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Fig. 7. — Top Panel: The measured planet-to-star flux ratios compared to the band- 
averaged ratios from atmospheric models that incorporate extra optical absorbers placed 
near the 0.01 bar level. Three models shown here in orange, blue, and purple, have absorber 
opacities Ke= 0.0, 0.1, and 0.0 cm^ g~\ and redistribution parameters P„ = 0.1, 0.3, and 
0.5, respectively. 

Bottom Panel: The measured flux ratios compared to the predicted ratios from the best-flt 
atmospheric model, with k=0.05 cm^ g~^ and P„ = 0.1, and the absorber placed near the 
0.1 bar level, deeper in the atmosphere than for the other models. 
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HJD 


Orbital Phase 


Observed Flux 


Flux Error 


2,454,846.686842 


0.4320984781 


0.9939594865 


0.0052779813 


2,454,846.687166 


0.4322515726 


1.0044958591 


0.0052779813 


2,454,846.687398 


0.4324662685 


1.0090936422 


0.0052779813 


2,454,846.687722 


0.4326120615 


0.9946241379 


0.0052779813 


2,454,846.687942 


0.4328347445 


1.0004521608 


0.0052779813 



Table 2: Final light curve, calculated by averaging the point-by-point average of the results 
from analyses A and B. Note: Table 2 is published in its entirety in the electronic edition 
of the Astrophysical Journal. A portion is shown here for guidance regarding its form and 
content. 
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Stellar Parameter Value Unit Ref. 

Stellar mass 1-0 1^0:22 ^0 3 

Stellar radius l.^blto^^ Rq 3 

Star Te// 5950±150 K 1 

log^ 4.25±0.30 (cgs) 1 

[M/H] -0.3±0.25 dex 1 

Right ascension (J2000) 06 48 19 1 

Declination (J2000) -03 06 08 1 

Measured Planet Parameter Value Unit Ref 

Planet mass Mp I-OTIqjs 3 

Planet radius Rp l-45l[!;?^ Rj 3 

Transit epoch Tq 2454159.452879±0.000068 HJD 2 

Orbital period 1.5089686^[!:S|]^ day 2 

Semi-major axis a 0.0259+[!:|J[J^J AU 3 

Orbital inclination i 85.66lQ;4g degrees 3 

Mid-eclipse phase 0.5022l[5:°[!P 4 

Orbital eccentricity e cos u; 0.0035+o:oo42 4 

2.2 /xm eclipse depth 0.336± 0.042% 4 

Blackbody model-derived Planet Parameter Value Unit Ref 

Bond albedo Ab ^.m^t^^o 4 

Re-radiation factor / 0.450;[|:^| 4 

Blackbody model temperature 2454l;i7o ^ 4 

Physical model-derived Planet Parameter Value Unit Ref 

Bond albedo Ab 0.075±0.015 4 

Energy redistribution factor P„ 0.1 4 

Absorber opacity 0.05 cm^g~^ 4 

Absorber depth in atmosphere 0.1 bar 4 

Temperature at absorber depth 2200 K 4 



Table 3: Star and planet parameters. References. (l) lBarge et al.ll2008l . (2) lBeanll2009l . (3) 



